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Abstract 

We present a measurement of the branching fraction for the doubly-Cabibbo-suppressed decay 

— ^ K^tt'^, using 281 pb^^ of data accumulated with the CLEO-c detector on the '0(3770) 
resonance. We find B {D+ K+7r°) = (2.28±0.36±0.15±0.08) x 10'"', where the first uncertainty 
is statistical, the second is systematic, and the last error is due to the uncertainty in the reference 
mode branching fraction. 
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The Cabibbo- favored hadronic decays of the c quark proceed through c sWy, Wy 
ud (Wy a virtual boson). The doubly-Cabibbo-suppressed decays proceed through c — >■ 
dWy, Wy — >• us, and are expected to be suppressed by a factor \{VcdVus) / (ycsVud)\'^ ~ 2.5 x 
10"'^. The doubly-Cabibbo-suppressed decay D° K~^7i~ was first observed in 1994 
and its branching fraction is now known to good precision (±2.8%, relative 0). Its ratio 
to the Cabibbo-favored decay K^n^ is measured to be (3.76 ± 0.09) xlO~^ 01, in 

qualitative agreement with the simple expectations. Very recently BaBar has observed 0] 
a second D — > Kn doubly-Cabibbo-suppressed decay — ^ K~^7i^ (charge-conjugate mode 
K^-K^ implied also, throughout). Here we report confirmation of BaBar's result, with 
slightly better accuracy. These measurements can provide insight into the decay mechanisms 
for D Ktt: the validity of SU(3), and the roles of the annihilation, exchange, and color- 
suppressed spectator diagrams relative to the color-favored spectator diagram 0, 0- A 
more extensive picture will be provided by the measurement of the remaining two D Ktt 
doubly-Cabibbo-suppressed decays, K^ti^ and K^n^ . 

For this measurement, we have used a 281 pb^^ sample of e+e~ colliding beam events, 
collected at a center-of-mass energy of 3770 MeV. The events were produced with the CESR-c 
storage ring and detected with the CLEO-c detector. The data sample contains about 0.8 x 
10® D^D~ events (our target sample), one million D^D events, five million e+e^ — > mm, dd, 
or ss continuum events, one million e"''e" t^t^ events, and one million e^e^ jip' 
radiative return events (sources of background), as well as Bhabha events, /i-pair events, 
and 77 events (useful for luminosity determination and resolution studies). 

The CLEO-c detector is a general purpose solenoidal detector which includes a tracking 
system for measuring momenta and specific ionization (dE/dx) of charged particles, a Ring 
Imaging Cherenkov detector (RICH) to aid in particle identification, and a Csl calorime- 
ter for detection of electromagnetic showers. The CLEO-c detector is described in detail 
elsewhere 0013. _ 

The ip{3770) resonance is below the kinematic threshold for DDn production, and so 
the events of interest, e'^e" — > ?/'(3770) DD, have D mesons with energy equal to the 
beam energy. Having picked the particles being considered to make up a Z) meson, following 
Mark HI l9| we define the two variables: 



A-E = ^Ei — -Ebeam, (l) 
i 

and 

^bc = ^El,^-\Y^P,\^, (2) 

where Ei, Pi are the energy and momentum of each D decay product. For a correct com- 
bination of particles, AE will be consistent with zero, and the beam-constrained mass Mbc 
will be consistent with the D mass. 

In addition to K~^7i^, we have studied the singly- Cabibbo-suppressed decay 

7r+7r°, as a higher-rate decay possessing kinematics similar to K^tt^, and 

the Cabibbo-favored decay — /('"vr+vr"'", as a high-rate, low-background mode used for 
normalization. We distinguish between and vr^ using information from the RICH and 
dE/dx information from the central drift chamber. We identify vr^'s via vr^ — »• 77, detecting 
the photons in the Csl calorimeter. We require that the calorimeter clusters have a measured 
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energy above 30 MeV, have a lateral distribution consistent with that from photons, and 
are not matched to any charged track. We require that the 77 invariant mass be within 3 
standard deviations of the 7r° mass. The 7r° mass resolution is 5.4 MeV (Gaussian width a) 
for both D+ K+ti'^ and D+ vt+tt". The AE resolution is 14 MeV for D+ K+ti^, 15 
MeV for D+ 7r+7r°, and 5.6 MeV for D+ K-n+n+. The Mbc resolution is 1.90 MeV 
for D+ K+7i^, 1.96 MeV for D+ ^ vr+vrO, and 1.35 MeV for D+ K-n+n+. 

We select candidate combinations that have AE between —40 MeV and +35 MeV for 
K~^TT^ and 7r"'"7r'^, and between —20 MeV and +20 MeV for K^-n-^n^. These requirements 
correspond to roughly 3 standard deviations. The asjTiimetric cut for K~^tt^ and vr+vr^ is due 
to a low-side tail on 7r° energies, and the wider window is due to poorer energy resolution. 
To study background, we select combinations with AE between —100 and —50 MeV, and 
between +45 and +100 MeV (+50 and +100 MeV for K~7r~^7r~^). When an event contains 
more than one K~^Tr^ combination that passes our AE requirement (a 1.4% occurrence), 
we choose the combination with AE value closest to zero. Multiple candidates per event 
for 7r"'"7r° and for K'n'^n'^ are comparable in frequency, and are removed by the same 
procedure. Thus, we allow only one candidate per event per decay mode per D charge. For 
those multiple candidate events that contain a real — >■ K~^7r^ decay, Monte Carlo studies 
indicate that our algorithm for picking the "best candidate" gets the right one 2/3 of the 
time. Because the algorithm uses AE only, and our procedure for extracting yield uses a fit 
to Mbc, the algorithm introduces no bias. 
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FIG. 1: Mbc distributions of D+ K-7r+7r+, D+ vr+vr^ and ET+vr". The points are 

obtained by selecting the AE signal region, the shaded histogram is from the AE sidebands, and 
the lines are the fit described in the text. 



The Mbc distributions for candidate combinations are shown in Fig. ^ The normalization 
mode — > K'n'^ir'^ is essentially background-free. The — »• tt+tt" mode background 
is well described by the distribution obtained from the AE sideband, as is that for the 
— s> K^Ti^ mode. There is a clear peak in K^tt^. 

Our Monte Carlo studies indicate that 80% of the background to — > /T+tt" comes 
from continuum events, 11% from DD events, 8% from radiative return events, and 1% 
from r-pair events. The AE requirement cleanly separates K'^n^ and tt+tt^ 

decays, so there is no cross-talk between these modes. There is no evidence for peaking 
backgrounds. 

We perform an unbinned maximum likelihood fit to extract signal yields from the Mbc 
distributions. For the signal, we use a Crystal Ball line shape [10], which is a Gaussian 
with a high-side tail. As Monte Carlo studies show that K^n^ and tt+tt^ 

have the same signal shapes, we have determined the line shape parameters (Gaussian peak 
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location, Gaussian width, point at which high-side tail begins) from the ti^ti^ Mbc 

distribution, and used them in the fit to the — >■ K^n^ Mbc distribution. We have varied 
the shape of the high-side tail as part of the systematic error study. For the background, 
we use an ARGUS function with shape parameter determined from the sideband 
Mbc distribution, high-end cutoff given by -Ebeam, and normalization determined from the 
fit to the AE signal region. Monte Carlo studies demonstrate that the shape parameter 
determined from the AE sideband correctly describes the shape of the background in the 
AE signal region. We have also performed a fit with the ARGUS shape parameter free in 
the fit, and obtained essentially the same result. 

TABLE I: The efficiencies (from Monte Carlo, but corrected for Tr'^-finding (see below)), fit yields 
from data, and branching fractions from data. Only statistical uncertainties are included. 



Mode 




Signal yield 


B (%) 


D+ K-i:+TT+ 


52.16 ± 0.16 


79612 ± 291 


9.51 (Input) 


D+ ^ vr+vrO 


47.65 ± 0.15 


964 ± 54 


0.1326 ± 0.0075 




42.30 ± 0.14 


148 ± 23 


0.0228 ± 0.0036 



Results of the fits are shown in TableHl Also given in Table|l]is the detection efficiency for 
each mode, and the branching fractions obtained for — >■ h^tt^ and — >■ K^n^. Those 
branching fractions are obtained by measuring the respective efficiency-corrected yields rel- 
ative to that for K^ii^n^, taking that branching fraction as (9.51 ± 0.34)%, which 
is taken from the 2006 Particle Data Group (PDG) value The branching fraction for 
—>■ 7r"'"7r*^ is in good agreement with our previously-published branching fraction using 
the same data set, (0.125 ± 0.006 ± 0.007 ± 0.004)% 12]. We emphasize that these resuhs 
are not independent, and the value in this paper should not be used in place of the previous 
result. 

We have considered many sources of systematic error to the —* K^ir^ branching frac- 
tion, including: signal Monte Carlo statistics, track-finding efficiency, vr^-finding efficiency, 
particle identification, the AE requirement, final state radiation, and the uncertainty from 
our fitting procedure (background shape, signal shape). The only ones greater than 1/10 of 
the statistical error are vr^-finding efficiency, background shape, and signal shape. 

The Monte Carlo simulation of the calorimeter response to photons is imperfect, partic- 
ularly in those angular regions where there is considerable material between the interaction 
point and the calorimeter. The Monte Carlo simulation overestimates the efficiency for de- 
tecting vr^'s. Various data-Monte Carlo comparisons suggest a correction factor of (0.95 ± 
0.04), which we apply. 

The background shape is determined by a fit to the AE sideband data. The error on the 
shape parameter thus obtained translates into a ±4.4% relative error in the — > K^t:^ 
branching fraction. The signal shape is determined by a fit to the tt+tt^ signal. 

Uncertainty comes from the determination of Gaussian width (cr), and the point at which 
non-Gaussian tail sets in (a). We have obtained the error ellipse in the determination of these 
two parameters, and noted the variation in fitted — >■ K^ir^ yield as one travels around 
this error ellipse. In this way, we obtain a relative systematic error of ±2.6%. Note that 
both the background shape uncertainty and signal shape uncertainty are really statistical 
errors, hence will decrease as additional data are taken. 
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We have also considered systematic errors to our normalizing mode, K~7t~^'k~^, 
i.e., to the yield and to the efficiency. Because this mode is essentially background-free, 
background shape and signal shape contribute negligible errors. Kaon particle identification 
tends to cancel in the ratio to K^tt^. Pion particle identification efficiency is well- 

modeled by Monte Carlo simulation. Track-finding efficiency - 3 tracks in normalizing mode 
vs. 1 track in signal mode, with 0.7% uncertainty per track - is the largest error, and is less 
than 1/10 the overall statistical error (1.4% vs. 16%). 

Our final result is 

B (d+ ir+7r°) = (2.28 ± 0.36 ± 0.15 ± 0.08) x 10'^, 

where the first error is statistical, the second error is systematic, and the third error is from 
the uncertainty in the K^n^n^ branching fraction, (9.51 ± 0.34)% 13], used as the 

normalizing mode. 

Our result is in good agreement with the only other measurement of this branching 
fraction, BaBar's recent B {D+ K+ti^) = (2.52 ± 0.47 ± 0.25 ± 0.08) x 10"^ |3|. It can be 
converted to a width, using the PDG value for the D+ lifetime ((1040 ± 7) xlO^^^ s) 0, 
and compared with the width for doubly-Cabibbo-suppressed D*^ decay D° — > K~^7r~ , using 
the PDG values for the D° K+tt" branching fraction ((1.43 ± 0.04) x 10"^) and 
lifetime ((410.1 ± 1.5) x 10"^^ s) 0. In this way we obtain 



r(DO ^ K+Ti-) B{D^ K+Ti-) X td+ ■ 

The spectator model diagram, expected to be the dominant contribution, predicts 1/2 for 
the ratio. Annihilation and exchange diagrams, which contribute differently to the two 
decays, can shift the ratio away from 1/2. Our result, and the BaBar result 0, suggest that 
such a shift is small. 

In summary, we have measured the branching fraction for K^tt^ to be (2.28 ± 

0.36 ± 0.15 ± 0.08) X 10~^, in agreement with the only other measurement of that branching 
fraction, and of comparable accuracy. 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. D. Cronin-Hennessy and A. Ryd thank the A. P. Sloan 
Foundation. This work was supported by the National Science Foundation, the U.S. De- 
partment of Energy, and the Natural Sciences and Engineering Research Council of Canada. 
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